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Abstract 

We study the dynamics of quenched fundamental matter in M = 2* supcrsymmetric large N c 
SU(N C ) Yang-Mills theory, extending our earlier work to finite temperature. We use probe 
D7-branes in the holographically dual thermalized generalization of the Af = 2* Pilch-Warner 
gravitational background found by Buchel and Liu. Such a system provides an opportunity 
to study how key features of the dynamics are affected by being in a non-conformal setting 
where there is an intrinsic scale, set here by the mass, m H , of a hypermultiplet. Such studies 
are motivated by connections to experimental studies of the quark-gluon plasma at RHIC and 
LHC, where the microscopic theory of the constituents, QCD, has a scale, Aqcd- We show that 
the binding energy of mesons in the N = 2* theory is increased in the presence of the scale mn, 
and that subsequently the meson-melting temperature is higher than for the conformal case. 



1 



1 Introduction 



In making further progress in strengthening the connection between the string theory based holo- 
graphic techniques and experimentally accessible strongly coupled systems such as the quark-gluon 
plasma studied in heavy ion collisions at facilities such as RHIC and LHC (for reviews see e.g., 
refs. [H [2] and references therein), better understanding is needed of systems well away from the 
conformal M = 4 theory. While it is remarkable that finite temperature N = 4 theory (at large 
N c , for SU(N C ) gauge group) already has several properties in common with the quark gluon 
plasma (e.g., the strikingly small shear viscosity to entropy density ratio for its hydrodynamic 
behaviour [31 Hj), the introduction of temperature, while indeed breaking conformal symmetry and 
supersymmetry, introduces only one scale into the problem, in terms of which all other thermal 
properties are determined. This is evident in the dual geometry, which is AdSs-Schwarzschild 
(times S 5 ) [3 [6], with an horizon at a radius that sets the temperature. That system (and its dual 
geometry — in local AdS coordinates) has no deconfining phase transition at some non-zero T c , 
and is in its "high" temperature phase for any T > 0. 

In contrast, the microscopic theory of nuclear matter, quantum chromodynamics (QCD), 
dynamically generates a natural scale, Aqcd, as a result of asymptotic freedom. Aqcd, ultimately 
helps determine the basic size of the bound states of quarks and gluons in the theory, setting a 
finite transition temperature for deconfinement to the quark-gluon plasma (QGP) phase. This 
extends to mesons in the theory, even after the QGP has formed. The properties of the spectrum 
of "quarkonium" , as a function of quark mass, and the subsequent melting of the mesons at tem- 
peratures T me it > T c , is a subject of considerable interest, both experimentally and theoretically, 
not the least because quarkonium suppression could be a valuable diagnostic probe in QGP studies, 
were it sufficiently computationally understood. 

While there are many active discussions of the extraction of T me it in the literature (see e.g., 
ref.jjTJ's review section and references therein) the analogous quantity in (so-called "top-down") 
holographic duals is easy to extract (at large N c , and in a quenched limit where the quarks do 
not back-react on the physics, so things are much simpler). There, meson melting is a first order 
phase transition [8, 9J, [10] from D7-brane probe [TT] fluctuations with boundary conditions leading 
to infinitely long-lived bound states (the mesons) to fluctuations of D7-branes that end on a black 
hole, where there the boundary condition gives states that decay by falling into the black hole (the 
"quasinormal" modes — the melted mesons). Nevertheless, the placing of D7-branes in the AdSs- 
Schwarzschild background (in the zero back-reaction limit) will still result in a melting temperature 
that recalls too much of the conformal nature of the parent theory. The /3-function of the resulting 
N = 2 theory at zero quark mass is again vanishing (since it is of order Nf/N c and here N c >> Nf 
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in the probe/quenched limit) and so T me i t will be determined only in terms of the quark masses, and 
since there is no intrinsic non-thermal scale like Aqcd 5 will be lower than what would be expected 
experimentally. 

To make progress it is highly desirable to move away from the neighbourhood of a conformal 
system and cleanly study meson melting in a holographic dual of a gauge theory with a natural 
scale present already at T = 0. The mesons size and binding should depend upon this scale (in 
addition to the bare quark masses) and so the gauge theory would have a T me i t that is higher than 
if it had no scale. 

This is the journey begun in our previous paper [12], and continued in the present one. The 
gauge theory in question is N = 2* , obtained by making massive an N = 2 hypermultiplet within 
the TV = 4 vector multiplet, and then flowing (already starting at strong coupling) to the infra-red 
(IR) to study the low energy physics. Now, there is no finite Aqcd in the IR for this theory, and 
instead the natural scale is set by the hypermultiplet mass mg, appearing in the dual geometry as 
a structure that is an example of the enhangon mechanism |13j, a beautiful combination of large N c , 
quantum effects, and strong coupling. The natural coupling, A = gy M N c , already starts out large 
for this theory (in order to be accessible as a smooth gravity dual at least asymptotically) and 
so dimensional transmutation generates no natural scale Aqcd in. the IR. This is the scale that 
normally generates (through quantum effects) the separation between monopole/dyon points on 
the Coulomb branch a la Seiberg— WittenflU 115] . At large N c , there generically would be N c such 
points and since Aqcd is vanishing here they are densely packed together, spreading out to form a 
circle, or in the natural coordinates, a line segment (see refs. |16L I17j ). The radius of the circle (or 
length of the line) is set by mn , so the hypermultiplet mass mn plays the role of Aqcd for us in the 
J\f = 2* theory, where the dual geometry is seen to be describing part of the Coulomb branch of the 
theory. In the dual geometry (known as the Pilch- Warner geometry[18j), the Coulomb branch is 
the place where the potential of a probe D3-brane vanishes, and the enhangon is the large N c locus 
where their tension drops to zero|161 117]. More generally, in the natural supergravity coordinates, 
the hypermultiplet mass ran sets a radius in the core of the ten-dimensional geometry where the 
metric becomes singular as one of the fields, x-> which in the ultra-violet (UV) sets the mass of the 
hypermultiplet, diverges. 

In our previous paper [12] we holographically studied quarks and mesons in the N = 2* 
theory at T = by exploring D7-brane embeddings in the dual ten-dimensional Pilch-Warner 
geometry. As we shall discuss later, our results already contained the seeds of the anticipated finite 
temperature result — a higher T me n than for the conformal theory — and our explicit study of finite 
temperature embeddings in this paper confirm this, also uncovering a great deal of interesting and 
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useful detail. We use the finite temperature geometry generalizing the Pilch- Warner geometry that 
was found in ref.|19j. and which has been studied extensively in the context of models of QGP-like 
dynamics in refs. |T9l EUJ [21] [22] [23]). Our explorations allow us to study a family of D7-brane 
embeddings and extract from them the behaviour of the meson melting temperature T me i t as we 
vary quark mass and ma- 
in section [2j we briefly review the zero temperature background to establish conventions 
and notation. In section |3| we review the thermal version of the N = 2* background. In section [4] 



we probe the background with D7-branes, extracting in subsection 4.1 the condensate of the A = 3 
operator given by [23] : 



(O) =i^ + q (m q + V^ 3 ) q ] + q ] (rn q + v^s) q + h.c. , (1) 

where 4>3 is the complex scalar in the N = 2 vector multiplet and the bare mass m q is the source 
of the operator. The fields (ft and (q, q) are the fermionic and bosonic components of the quark 



multiplet. We discuss meson spectrum and the meson-melting temperature in subsection 4.2, and 
conclude in section [5] 

2 Zero Temperature J\f = 2* and Dual Extremal Geometry 

The matter content of the J\f = 4 supersymmetric Yang-Mills theory consists of a gauge multiplet 
containing the bosonic fields (A^, Xj), j = 1, . . . , 6, where Xj are real scalars transforming as the 
6 of SO(Q), and the fermionic fields (Xj), j = 1, . . . , 4, which transform as the 4 of SU (4). Writing 
this matter content in terms of N = 1 superfields, the theory has one vector supermultiplet (^4^, A4) 
and three chiral multiplets &j = (Xj, 4>j = ^2j-i + 1X2 j), j = 1, . . . , 3. We can make two of the 
chiral multiplets massive (1^1 and $2, say, with equal mass), preserving only an M = 2. 

In the dual gravity picture, this will correspond to turning on two real scalars a and \ with 
conformal dimension A a = 2 and A x = 3 respectively. The theory described by the resulting flow 
to the infrared (IR) is called the M = 2* theory. The flow can be holographically described in 
terms of five dimensional gauged supergravity, with action given by: 

S 5d = J d 5 xV^G (ll - 12 (da) 2 - 4 (d X f - 47>) , (2) 

where the potential V for the scalars x an d a is: 



V = 9 - 
16 



1 fdW\ 2 fdW^ 2 
+ 



3 V da J \ dx 



a 2 1 

^-W 2 , with W = -e~ 2a - -e 4a cosh(2x) . (3) 



The constant g is related to the AdS radius R, g 2 = 4/B 2 . 
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In the zero temperature (supersymmetric) case, the problem of solving for x( r )> p( r ) an d 
A(r) where the latter is given by the (Einstein frame) metric ansatz: 

dsl = e 2A {-dt 2 + dx 2 ) + dr 2 , (4) 

reduces to solving the following first order equations |18j : 

dp 9 (l P 4 c«W) , (5) 



dr 3R \p 2 

dA 2 ( 1 1 4 . A 

Partial solutions are given by |18j : 

ko 2 

&A = ■ S ' P S = cosh ( 2 x) + sinh 2 (2 X ) (In (tanh( X )) + 7) • (8) 
smn(2xj 

Here k = mj{R, where mjj is the mass of the hypermultiplet and 7 is a constant corresponding to 
different slices through the moduli space of M = 2* in the IR. 

Note that we can have a complete analytical solution for the metric in terms of the field x 
by making the following transformation: 

dx = -fdr ->■ dr 2 = . d X 2 , 9 

dr /9 8 smlr(2x) 

such that the AdS boundary is at x = 0- in t ne interior % — > 00 , and the apparently singular 
physics there was made sense of in terms of the Coulomb branch of the gauge theory using D3- 
brane probe techniques \T7\ PT6]. (For example, the enhangon locus appears there, as mentioned in 
the introduction.) 

Beyond that, we can solve for x numerically, and it is convenient to do this in a coordinate 
system given by: 

z= Z - = e~*l R . (10) 

It was shown in ref. [19] that in terms of this coordinate, the fields have an expansion near the AdS 
boundary given by (for 7 = 0): 



X = k'z 



,2,2 1 1 , 4, , lA ,J 7 . 10 20 



1 + k 2 z 2 [- + - ln(kz) ) + k^z^ ( + — ln(M) + — ln(kz) 2 ) + 0(k b z b ]n(kzy) 

90 3 9 



p = 1 + k 2 z 2 Q + I ln(H)^ + k 4 z 4 (J^ + 2\n(kz) + | ln(A-- ) 2 ) — 0{J/'-J' \n( l:z/' ) . 

A = - ln(2l) - -k 2 z 2 - k 4 z 4 (- + — ln(kz) + - ln(A;z) 2 ) + 0(k 6 z 6 In(kzf) . (11) 
3 V 9 9 9 ' 
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3 Non— extremal Background 



We follow the analysis of ref. |22j . which we repeat here for completeness. Consider this form for 
the metric: 

dsl = - Cl (r) 2 dt 2 + c 2 (r) 2 dx 2 + dr 2 . (12) 
With this ansatz the equations of motion for the scalars give, in terms of the scalar potential P: 

a"(r) + d r (In ( Cl (r)c 2 (r) 3 )) a'(r) - = , 

X "(r) + d r (In ( Cl (r)c 2 (r) 3 )) X '(r) -l^-=0, (13) 

and the Einstein equations give: 

c'((r) + d r (31nc 2 ) c'^r) + ~ Cl (r)P = , 
4(r) + d r (In ( Cl c|)) ci(r) + ^c 2 (r)P = , 
and a (r) 2 + ^x'(r) 2 -^V-^d r (In c 2 (r))d r (In Cl (r)c 2 (r))=0 , (14) 

where for the last equation we used the previous two equations to simplify its form. 

Since we will be solving these equations numerically, the AdS radial coordinate r turns out 
not to be convenient, since it runs to infinity (the location of the AdS boundary). Therefore, we 
define a coordinate x(r) satisfying: 

1 - *(r) = — 7—r , (15) 
c 2 (r) 

such that x G [0, 1], with x = being the AdS boundary and x = 1 being the event horizon. An 
important thing to note is that this change of coordinates will only work for the non-extremal case, 



since for the extremal case, we have c\(r) = c 2 (r) and so x = for all r. Using equation (14), we 
can derive an expression for dx/dr entirely in terms of x: 

dx =± I 2Vc 2 (x) 2 (l-x) ~ 

dr y (l-x)c 2 (x) 2 (6a'(x) 2 + 2 X '(x) 2 )-3c/ 2 (x)(-c 2 (x) + 2(l-x)d 2 (x)) ' 1 ' 

In this new coordinate system, the equations of motion reduce to: 

4(x) - 5^ + + \ c 2 ( x) ( 3a ' (x) 2 + X '(x) 2 ) = , 

C 2 {X) 1 — X 6 

a'(x) 1 OP 
a {x) -—x-6 Gxx da- = > 

and X"(x)-^1- 1 -G xx ^ = 0, (17) 
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where 

(1 - x) c 2 (x) 2 (6a'(x) 2 + 2 X '(x) 2 ) - 3c' 2 (x) (-c 2 (x) + 2(1 - x)c' 2 (x)) 

2Pc 2 (x) 2 (1 - x ) • [L *> 

Note that in order to derive these results, one needs the quantity: 

4c' 2 (x) (-2(1 - x)c 2 (x) 2 (3a'(x) 2 + X '(x) 2 ) - 3c 2 (x)c' 2 (x) + 6(1 - x)c' 2 (x) 2 ) 

P(l-x)c 2 (x) 3 ' [ ' 

where we have used the equations of motion to simplify the expression in terms of only first 
derivatives of the fields. It is worth recalling here the connection to the fields of the previous 
section: 

p{x) = e a W j C2 ( x ) = eA(x) _ (20) 

3.1 Non— extremal AdSs 

Let us consider some special cases to orient ourselves. Consider the case of a = x = 0. The most 
general solution for A(x) is: 

A(x) = lnai- -1n(x(2 - x) + a 2 ) ■ (21) 

We can take a 2 = 0, and the metric is now given bjj^j 

dsl = a\ (x(2 - x)r 1/2 (- (1 - x) 2 dt 2 + dx 2 ) + ^ (2x - x 2 )" 2 dx 2 . (22) 

We can return this to the perhaps more familiar form of AdSs Schwarzschild by the following 
coordinate transformation: 

x = l-<f-^, (23) 
and the choice of a\ = b/R = (irTR) such that we have: 



3.2 Non— extremal Pilch— Warner Geometry 

Motivated by our result above for non-extremal AdSs, we consider an ansatz for the fields of the 
following form: 

A(x) = ln(«) - - ln(x(2 - x)) + a(x) . (25) 



1 There appears to be a typo in the g xx term for this equation in ref. [22] . 
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Near the AdS boundary, the fields have the following leading behavior: 



p(x) 
X{x) 

a(x) 



1 + x 1/2 (pio + pu ln(x)) + 0{x) , 

Xo* 1/4 fl + x 1 ' 2 ( X w + \xl H*)) + O(x)) , 



Xlx^ + Oix) . 



(26) 



The parameters (pi 1, Xo) are related to the bosonic and fermionic masses (mg and mp respectively) 
of the components of the (would-be) N = 2 hypermultiplet) [22j : 

. 2 24vr 2 



V T 



V2 



(27) 



where we have defined = a(l). Furthermore, the parameter k is related to the temperature via 
1221: 



-30^ 



i?7T 



(28) 



Therefore, to proceed, we fix the values of (k, /Oi,i,Xo) an d search the values of (/9i,o,Xl,o) that 
give us regular solutions (ie. solutions for (p(x), x( x )> a ( x )) that are have vanishing first derivative 
at the event horizon). We focus our attention on a particular type of solution, which we label as 
ref. [22] does as a "supersymmetric deformation" with mp = mg = mu, which corresponds to 
taking: 

Xo = 6/5i,i • (29) 

We present our numerical solution for this case in figure [TJ and they appear to agree very well with 
the results of ref. |22j . We will not consider cases of the type mp 7^ ms in the rest of this paper. 

Pl,0 



Xo 




7).05 0.10 0.15 
(a) SUSY Deformations 



0.2i 



(f 




(b) SUSY Deformations 
Figure 1: Non-Extremal Pilch- Warner solutions. 
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4 Quarks and Mesons from D7— branes 



The five dimensional physics may be oxidized to a ten dimensional system (more suitable for 
discussion of the D-brane probing) within the type IIB supergravity [18] , with the ten dimensional 
metric (in Einstein frame) given by: 



(30) 



dsj = tfdsl + ^ (V W + p 6 cos 2 9 (JL + °*±&) + sin 2 e^j , 
where: 

Xi (r, 9) = cos 2 9 + p(rf cosh(2x(r)) sin 2 9 , 
X 2 {r,9) = cosh(2x(r)) cos 2 9 + p{rf sin 2 9 , 

ft 2 = Mi^, a 2 = 4 = 2jR 2 ( 31 ) 

c = cosh(2x) , p = e a , 
and there is a deformed S 3 with 5^7 (2) invariant 1-forms 



| (da + cos ifid/3) , 

2 (— sin adijj + cos a sin^d/?) , (32) 



(71 : 

and o"3 = | (cos ad?/> + sin a sin i/)d/3) 



and the dilaton is given by: 

e = ^ — . (33) 

cX\ sin (f> + X2 cos 2 cj) 

There are several other fields, but we will not need them for the embeddings we will choose, following 
our work at T = in ref.[12], to which we refer the reader for details. Following our analysis in 
ref. [12], the relevant action for us is: 



S D7 = -fPr J A eV-det (P[G] ab ) , (34) 
with a D7-brane embedding given by: 

£ a = x a ,a = 0, ...,7 , 9 = 9(x), <p = vr/2, 3vr/2 . (35) 
Near the AdS boundary (x — > 0), the field 9(x) has asymptotic behavior given by: 

9(x) = x 1 ^ + x 3 / 4 (e 2 - ^x 2 lux) + . . . . (36) 



To extract the bare quark mass and the condensate, it is convenient to express equation (36) in 
terms of the coordinate z of equation (10). In order to do this, we compare the expansion of e A 



near the AdS background for the zero temperature background (which is in terms of z) with that 
of the finite temperature background (which is in terms of x), since near the AdS boundary, the 
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asymptotic behavior of both backgrounds should be the same. Using these two expansions and 
setting them equal, we can solve (iteratively) for x(z), 



x (z) = 8k 4 z 4 + 



128^ 
9~ 



Xo K z + • • • • 



Plugging the expression for x(z) into the asymptotic equation of 9(x), we get: 



^2 + ^0X0-^0X0 ln ( — 

» 3 vxo 



9(x(z)) = z—9 + z 3 ( — 
Xo VXo 

If we make the identification: 

9 = -flb , 02 = f-V (fc + jjfloag - Jfloxgln f- 
Xo VXo/ V 9 3 VXo 

and use the expression for the bare quark mass and the condensate |12j : 

m q = ^§ , (W) = 2ir 2 R 5 fi 7 (-29 2 + ^ - h 2 9 ln^ ) = 2it 2 R^ 7 C 
we find that the bare quark mass and condensate are given by: 



m„ 



2vra' VXo/ 
where we have used that k/xo = 2 3 / 4 k. 

4.1 Solving for the Condensate 



~ 29 ° + l e3 ° ~ t 60 ** ~ t 9ox2 ° Heo) I = h "'' c ' 



(37) 



(38) 



(39) 



(40) 



(41) 



We show the results for our condensate in figure 2(a) . We find that there is a divergence for large 



quark mass, a situation that we have encountered in the zero temperature case |12j . In the zero 
temperature case, the divergence can be attributed to the numerical artifacts arising from solving 
for 9 on top of a numerical solution for \- We can perform a similar subtraction as in the zero 
temperature case to get the result in figure 2(b)| We see that there is still a divergence that is 
unaccounted for. This is unsurprising since for finite temperature, we have more numerics in the 
background (p and a). In particular, the asymptotic behavior of a(x) is what gives the coefficient 



of both the linear term and linear-log term in equation (41). Therefore, we should expect that 



numerical artifacts would appear as exactly linear and linear-log terms. This is indeed what we 



see in figure 2(b) . We can subtract off this divergence by fitting the asymptotic behavior, and get a 
physical result as shown in figure 2(c) We give plots in terms of physical quantities in figure [2(d)] 
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c 




(a) Raw Result (b) Result after T = subtraction 



< ifr ifr > 




(c) Physical Result (d) Several Condensate Curves 



Figure 2: An example of the results (for 1 ^ — 0.3193) and subsequent subtraction scheme results. In (a), 
we have our raw results. In (b), the solid line is a curve of the form ax + bx In x that is fit to the condensate 
for large bare quark. After subtracting this curve, we end up with the result in (c). In (d), we give curves 
for ' T £f- = 0,5.12418,6.93171 in ascending order. Throughout, the black hole embeddings extend from the 
left, starting at the origin, while the Minkowski embeddings extend to the right. They meet at a kink that 
hides a region of multivaluedness discussed later in the text. 



4.2 Meson Melting Temperature 

We now wish to study the effect of the hypermultiplet mass on the melting temperature of the 
mesons, as discussed in the introduction. The meson is associated with 7-7 strings and its mass 
can be determined by considering fluctuations of the D7-brane embedding [25J. Recall that the 
Minkowski embeddings are related to the phase where we have mesons (quark bound states) in 
our gluon plasma, whereas the black hole embeddings are related to the phase where the mesons 
have melted in the plasma (the finite lifetime quasinormal modes arising from in-falling boundary 
conditions at the horizon). We may predict what might happen based on revisiting the zero 
temperature result. We revisit our results from ref. |12j . but present them in a different fashion in 
figure [3} 

Recall that the meson mass should in general be given by: 

M m = 2m c - E b , (42) 
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(a) 



2 4 6 8 

Meson mass for fixed bare quark mass. 




1 2 3 

(c) Detail of part (a). 




2 4 6 8 

(b) Meson mass for fixed constituent mass. 




2na i 



1 2 3 4 B? m„ 

(d) Detail of part (b). 



Figure 3: Meson spectrum in the extremal Af = 2* background. The red, dashed line is the meson spectrum 
in the Af = 4 background, whereas the blue, solid line is the meson spectrum in the M = 2* background. 



where m c is the constituent quark mass and Ef, is the (positive) binding energy. The constituent 
mass is calculated from the Nambu-Goto action per unit time for the string that hangs (in the x 
direction) from the "end" of the D7-brane at 9{x = xq) = n/2 to the event horizon at x = 1: 



m c 



—— [ dxe^ ,2 y/-g tt gxx 
2vra' J xn 



(43) 

In these scenarios, the meson mass is a factor of A -1 / 2 relative to the constituent mass, so the meson 
is said to be deeply bound [26] . For this reason, we cannot directly extract the binding energy for 
the meson at a given constituent mass, however we can make relative comparisons between the 
constituent mass in Af = 2* versus the constituent mass in the M = 4 for fixed meson mass. Our 
result in figure [3] suggests that the same mass meson in = 2* has a higher constituent mass than 
the meson in Af = 4. This in turn suggests that the (positive) binding energy in N = 2* is larger. 
Therefore, we would predict that in the thermal case, we should expect a raising in the melting 
temperature. 

To proceed, we can study the meson mass associated with different Minkowski embeddings 
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as follows. We consider fluctuations of the field 0{x) such that: 

9(x) = 6 {x) + 2vra'$(t, x) , (44) 
where 6q(x) corresponds to the embedding solution. We can expand the Dirac-Born-Infeld (DBI) 



action given previously (see equation (34)) to quadratic order in $(i, x) and derive the equations 
of motion for x) from the resulting Lagrangian. We take as an ansatz for the time dependence 
of <&(i, x): 

<5>(t,x) = e-' l Ti M 0{x) , (45) 

such that the meson mass (associated to this fluctuation of the D7-brane) is given by: 

2k 

M = —Co . (46) 
R 

Note that the factor of 2 is because of the factor of 1/4 associated with Gtt m the UV. The allowed 
values for ui are computed by requiring that 4>(x), which has expansion in the UV given by: 

4>{x) -> x 1/4 0i + ^ 3/4 ^2 , (47) 

only has normalizable modes, i.e. <j)\ = 0. We present some solutions for the meson mass in 
figures [5j Several important features include the crossing over of the meson mass ratio M m /m q 
below that of thermal J\f = 4 as we approach T/m q = in figure [5(a)j This was already observed 
in the zero temperature results in figure 3(a)| 



The meson melting transition occurs near where the meson mass is zero (the far right of 
figures [5]), but not exactly. To find the exact location (in terms of the quark mass), we need 
to calculate the free energy of the system and find where the jump (determined by minimizing 
the free energy) between the black hole and Minkowski embeddings occurs, inside the region of 



multivaluedness of the condensate vs. mass curves of figure 2(d) We have enlarged an example of 
such a region in figure [4j (See e.g., ref. [S] for examples of this extraction of the temperature.) 

We do not have numerical control of the free energy. However, where the meson mass drops 
to zero (corresponding to the end of the Minkowski embeddings) provides a good estimate of the 
melting temperature especially in terms of the bare quark mass. This is because the crossover 
region between the black hole and Minkowski embeddings is a very small part of the overall curve 
(this becomes even more accurate as ran increases). In terms of the constituent mass, we believe 
that the general behavior is correct although we should not take the exact numerical values here 
too seriously since the constituent mass rises very rapidly near the zero mass meson. With this 
in mind, we can systematically find an estimate of the melting temperature, and we present our 
results in figure [6j As predicted by the T = results, the melting temperature increases as we 
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<lf/lp> 

2 if R 5 fi 7 
0.35 r ■ 



0.30 
0.25 
0.20 
0.15 



2a' m q 



1.620 1.625 1.630 1.635 1.640 1.645 1.65Q? 2 T 

Figure 4: An enlarged portion of a condensate vs. mass curve showing the multivaluedness in the region 
where the melting transition occurs. The black hole cmbcddings enter from the left, while the Minkowski 
embeddings enter from the right. 



increase mn- Notice that the melting temperature saturates for very large run, a fact that we 
discuss further in the conclusions. 




R T 



0.6 2 a m 



(a) Meson spectrum for fixed bare quark 
mass 




(b) Meson spectrum for fixed constituent 
mass 



Figure 5: Solutions for the meson spectrum for different values of mn/T. The curves correspond to the 
following values of m H /T: (Red dashed, 0), (Green dot-dashed, 3.74), (Blue solid, 5.12). 



5 Conclusions 

In a very clean setting, we have been able to demostrate how the move to a non-conformal system 
changes the quark and meson dynamics in an important way. The instrinsic scale, mn, of the 
N = 2* theory affects the temperature at which mesons of a given mass melt (an important 
phenomenon in e.g., the experimentally accessible quark-gluon plasmas studied at RHIC and LHC). 
This is because the quark bound states have their properties determined by not just the scale set 
by the bare quark mass, but also by m#. There is a strong analogy to QCD, with its dynamically 
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1 2 3 4 5 6 7 2a m q 

(a) Melting temperature for fixed bare quark 
mass 



R 2 T mch 
2 a m c 




500 1000 1500 2a m c 

(b) Melting temperature for fixed constituent 
mass 



Figure 6: The melting temperature for fixed bare and constituent mass. 



generated scale Aqcd, and indeed this was in part the motivation for our study. 

Of course, the M = 2* theory is still quite different from QCD, and so we should expect 
differences in the meson dynamics. For example, the scale of J\f = 2* is the mass of a field in 
the theory, while in QCD the scale is generated by dimensional transmutation. Indeed, we see 



from figure |6(b) that, for fixed constituent quark mass, T me it eventually saturates for high m#, 
a non-trivial dependence that one would not expect in QCD. We expect that this is due to the 
fact the the binding energy and constituent quark mass get their contributions from the massive 
hypermultiplet in different ways. The contribution to binding energy evidently becomes more 
suppressed at higher mn- If this were a weakly coupled scenario, an examination of the structure 
of the Feynman diagrams contributing to each quantity at a given order (keeping only planar 
diagrams, since we are at large N c ) might make this more manifest, but it is not clear if the physics 
of interest is at all visible in perturbation theory, since we are dealing with bound states. It would 
be interesting to investigate this further. 

Another key difference is that while we did have a scale, we were not working in a confining 
theory. The geometry tells us about the Coulomb branch of the J\f = 2* theory. To get access 
to quark and meson dynamics (including meson melting) in a confining theory, a study of how to 
holographically lift the N = 2* theory's Coulomb branch and trigger monopole condensation would 
be one of a number of interesting directions to pursue. 
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